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Abstract 

The quasiparticle density of states in a two-dimensional d-wave supercon- 
ductor depends on the orientation of the in-plane external magnetic field H. 
This is because, in the region of the gap nodes, the Doppler shift due to the 
circulating supercurrents around a vortex core depend on the direction of H. 
For a tetragonal system the induced pattern is four- fold symmetric and, at 
zero energy, the density of states exhibits minima along the node directions. 
But YBa2C306.95 is orthorhombic because of the chains and the pattern be- 
comes two-fold symmetric with the position of the minima occuring when H 
is oriented along the Fermi velocity at a node on the Fermi surface. The effect 

of impurity scattering in the Born and unitary limit are discussed. 
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I. INTRODUCTION 



It has become increasingly clear that including the circulating supercurrents outside the 
vortex core in a (i-wave superconductor can lead to a qualitative understanding - in the low 
magnetic field limit, near Hd - of the effect of a magnetic field H on various properties of 



of the data for H oriented perpendicular to the Cu02 planes including impurity scattering 



symmetry in the electronic thermal conductivity. Very recently results have been obtained 
by Vekhter et a/.llllli for the effect of an in-plane magnetic field on the quasiparticle density 
of states (QDOS) in a tetragonal two dimensional (i-wave superconductor. The present 
work extends this previous study to the orthorhombic (i-wave case and also includes in the 
calculations the effect of impurity scattering in both, Born's and resonant scattering limits. 

While most high oxides are tetragonal, YBa2Cu306.95 (YBCO), which is perhaps the 
most extensively studied of all the oxides, is orthorhombic because its crystal structure 
includes CuO chains oriented along the b well as two Cu02 planes per unit cell. 

Optical measurements reveal that the carrier density in the chains in optimally doped YBCO 
is of the same order of magnitude as it is in the planes and the ratio of the associated b- 
to a-direction plasma frequency is of the order of 2.0 Similar large in-plane anisotropies 
are observed in the infrared optical conduct ivity,lli d.c. resist ivity,0lll zero temperature 
penetration depth,liiJii and thermal conductivity.i^'0 

Various approaches have been taken to treat band anisotropy including a single tight 
binding band with different nearest neighbor hopping in a- and 6-direction, iitwoband&13 
as well as three band^lil models which include the single CuO chain layer as well as two 
Cu02 plane layers per unit cell with transverse hopping between them. All these models 
give the same qualitative features and describe reasonably well the a-b band anisotropy 
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of YBCO. Additionally, in an orthorhombic structure the superconducting gap will be a 
mixture0 of dominant (i^2__y2 symmetry as well as subdominant Sq and Sx^+y^ symmetric 
functions because they all belong to the same irreducible representation of the crystal 
point group. 

The simplest approach which contains the physics of the more complex calculations is a 
single infinite band with distinct effective masses along the a (ma) and b (mij) directions&B 
with a {d + s)-symmetric gap function on the cylindrical Fermi surfacciiS For simplicity 
we will adopt throughout this paper this approach and calculate the effect of an in-plane 
magnetic field H on the QDOS assuming a distribution of non overlapping vortex cores 
and account only for the circulating supercurrents around the vortices. We will also treat 
the effect of impurities in both, Born's and resonant scattering limits. Impurities were not 
considered in the previous calculations in tetragonal symmetry by Vekhter et a/.0 which 
were confined to the pure limit. Thus, our results with impurities in the tetragonal case are 
the first such results. 



II. FORMALISM 

We begin with a single free electron band in two dimensions with an ellipsoidal Fermi 
surface and the single particle energy in the normal state given by the dispersion relation 

{h = l) 

e^ = ^ + ^-eF (1) 
Irrix 2my 

where k = {kx,ky) is the momentum vector and ep the Fermi energy. In Eq. (|l|) rrix and 
rriy are effective masses with > rriy, i.e. chains are along the ky-axis which increases the 
electronic transport in that direction. A parameter 

a = "I^l—^ (2) 

nine + f^y 

is introduced and this is the single parameter which characterizes the band structure 
anisotropy in our simple model. Consistent with the above assumption for the energy band 
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dispersion relation is a mixed {d + s) -symmetric gap on the Fermi surface (an ellipse) of the 
form 



Ak = Ao(cos2(^ + s), (3) 

with Aq the gap amplitude. In Eq. is a polar angle measured from the kx-a.xis in 

the two dimensional copper oxide Brillouin zone giving the angle of k (momentum on the 
Fermi surface) as shown in Fig. |l[ The first term in Eq. (|^) gives the dominant d-wave 
part of the gap which, on its own, would have nodes on the main diagonal at (f = ±7r/4. 
The number s is a measure of the subdominant s-wave part. Both parts are allowed to 
exist in an orthorhombic system such as YBCO from group theoretical considerations. The 



s-component moves the nodes from Lp = ±7r/4 to if = tan^^ y (1 + s)/(l — s) off the main 
diagonal of the Brillouin zone. For positive values of s the critical ip defining the node is 
larger than ir/A while the opposite holds for negative values of s. 

In order to work with Eqs. (|l|) and (|^) it is convenient to transform to new coordinates 
P = {Px,Py) in which the Fermi surface is a circle (see Fig. |l]).0 The required transformation 
is given by 



Pi = kid — , i = x,y (4) 



which leads to (Fig. 



/ TTl 

tan0=J — -tamp, (5) 



where is now the angle of the momentum vector p in the transformed frame. Clearly, the 
electron dispersion relation (^ reduces to 

= Ph+Pl,y ^ P^ (g) 
2m 2m ^ ' 

with ffi = {rrix + rny)/2 the average band mass. Applying the same transformation to the 
gap (^ gives the order parameter in the p-frame&B 

A / a + cos 20 \ ^ „ , 
- \ 1 + a cos 20 I 



We show in Fig. ^ the gap Ap in the p- frame for several values of a (the band anisotropy) 
and s (the gap anisotropy). Frame (a) is for a = and s = and is the well known pure d- 
wave case for which plus and minus lobes have equal absolute magnitude and the zeros are on 
the diagonal. In frame (b) a = 0.4 and s = 0, i.e. the gap in the k-frame (laboratory frame) 
is of pure d-wave symmetry but the band structure orthorhombicity leads to a transformed 
gap with larger positive than negative lobes and the nodes are shifted upwards from ±7r/4 
and similarly for the other symmetry related angles. This trend is even more pronounced 
in frame (c) where a = s = 0.4. In this case a and s add constructively to further increase 
the positive lobes as compared to the negative ones, and the nodes are shifted to yet larger 
angles. Finally, frame (d) applies to a = 0.4 and s = —0.4. For this set of parameters the 
nodes in the p-frame remain on the diagonal but at the same time the negative lobes are 
dominant over the positive ones. We conclude that a negative value for s can in some sense 
counteract some of the effects of a positive a, but there is no complete cancelation between 
these two parameters which play quite distinctive roles. The nodes can be made to remain 
unshifted but the size of the lobes is changed over their a = s = values. 

In the laboratory frame (k, before transformation) the thermodynamic Green's function 
takes on the form 

X (i^n - Vi.(k)q,) To + + AkTi 

G(k,a;„;rj- — — — — , (Sj 

with ffk and Ak given by Eqs. (|l|) and (H) respectively. In Eq. (|) the r's are Pauli 2x2 
matrices and iun is the n'th Matsubara frequency, Un = {2n + l)vrT, n = 0, ±1, ±2, . . ., and 
T is the temperature. The Doppler shift caused by the presence of circulating supercurrents 
outside the vortex core (the region which we assume to dominate the physics) is given by 
Vir(k)qs, where Vi?(k) is the electron Fermi velocity on the ellipsoidal Fermi surface in 
k-space (see Fig. |l]) and which is given by 

v,.(k) = Vkek = ^x+^y, (9) 

Wiy 

with X and y unity vectors in the direction of the kx- and /c^-axis. Applying the transfor- 
mation @) to get Vi7'(k) in the p-frame leads to 



vf(p) = -r^ + — . (10) 



The momentum of the superfluid currents outside the vortex core (placed at the position 
r = in the laboratory frame) is assumed to be inversely proportional to the distance |r| = r 
from the core. For a magnetic field H in the Cu02 plane making an angle 7 with the kx- 
axis the supercurrents are in a plane perpendicular to H and therefore to the k^-ky plane. 
Assuming circular supercurrents, an approximation discussed in the paper by Vekhter et 
a/.,0 and denoting the vortex winding angle by /3 one obtains for the x and y component of 
the superfluid momentum 

(!s,x = -|qs| sin /5 sin 7, 

= -|qs| sin/5cos7, (11) 

so that 



VF(p)qs = ^^Iqsl sin/3vT+~a ( \ — sin0cos7 — cos0sin7 I . (12) 
rux y\l 1-a J 

1 / </>() 



If one measures the distance r in units of the intervortex distance R = a \/ where (po 



is the flux quantum and a a constant of the order of unity, we obtain 



VF(p)qs = — — sin Py/l + a i \/— ^— — sin cos 7 — cos sin 7 I , (13) 

p 1 - " 

with the normalized distance p = r/R and the magnetic energy 



with vp a suitably defined effective Fermi velocity^ which can also be made to contain the 



additional band factor a/1 + a of Eq. ([T3| ) and z/ is a parameter which measures the magnetic 
energy Eh in units of the gap. Consequently, we can write 

VF(p)q. = —H(j), 7) sin/3, (15) 
P 

with 
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7) = yY~~ — sin cos 7 — cos sin 7. (16) 

The QDOS at zero frequency obtained from the Green's function which includes the 
circulating supercurrents about the vortex core and which must also be averaged over the 
vortex unit cell is given by 

By the square root is meant that branch of the complex function which corresponds to a 
positive imaginary part. Eq. ( p^ can be reduced to an one dimensional integral form 



One good approximation to this integral is obtained by noting that the integrand in Eq. (|1 
is dominated by the nodal region and, thus, we get approximately 

nodes yrnj 

where the 0„ are the nodes of the transformed gap and f'{(f)n) is the derivative of this 
gap function at its nodes. These can be found at 



(19) 



tan0„ = ±, 



1 + «)(1 + .) 



-\ (l-a) {l-sY 

while in the laboratory frame, as we have noted before, the gap has its nodes at tanyj^ 



±y (1 + s)/(l — s). Furthermore, the derivative f'{(t)n) at the gap nodes can be worked out 
to be 



/'(0„) = 2(1 + s)^{l-s^){l-a^), (21) 

and hence an approximation for A^(0,7) follows from Eqs. (0), (|2l|) and the definition of 
/i(0, 7) given in Eq. (|l^) as 



iV(0, 7) = - E (22) 



TT 



nodes (l + s)J(l-s2)(l-a2) 



with distinct contributions from only two of the four nodes. 
For one set of two nodes we get 



COS( 



sin( 



ll{l-a){l-s) 



1 + as 



as 



(23) 



while for the other ones 
expression 

71 



IT — (pn- We therefore arrive at the approximate analytic 



1 1 + a 



;i-s2)(l-a2)l + sV 1 



a 



1 + a)\l cos 7 — (1 ~ ol)\\ - — ; sm7 



1 + as 



1 + as 



+ 



1 + COS 7 + (1 — ajW sm7 



1 + as 
1 



1 + as 



1 1 + a 



J(l-s2)(l-a2) 1 + sV 1- 



a 



xmax 



, / 1 + s 
1 + ajW cos 7 



'1 — q)\\ sin 7 



(24) 



where we have restored factors previously absorbed in the magnetic energy Eu that appeared 
explicitely in Eq. (^). In the tetragonal, pure d-wave case (a = s = 0) expression (| 
reduces properly to the known result^! 



A^(0, 7) = ^ max (| sin7|, | COS7I), 



TT 



(25) 



in which case Eq. 



is known to be an excellent approximation to the one dimensional 



integral of Eq. (|T8]) . 

When band anisotropy is included through a finite value of a or a subdominant s-wave 
part, characterized by the parameter s, the analytic expression (|2^) is not quite as good an 
approximation as ( p5l) was but it is still acceptable as can be seen in Fig. In this figure we 
compare with our exact numerical results for A^(0, 7) based on the one dimensional integral of 
Eq. (|TH). Here only band anisotropy is included with an effective mass anisotropy parameter 
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a = 0.2. The dashed hne is the exact result (|T8|), the sohd hne is based on Eq. (p^. The 
agreement is good. Note the two-fold symmetry in the pattern obtained, although, since 
s = 0, the gap in the laboratory frame has pure rf-wave symmetry with nodes on the main 
diagonal. It is the band anisotropy that breaks the four-fold symmetry of the tetragonal case. 
However, it is very important to realize that because of a 7^ 0, i.e. band orthorhombicity, 
the sharp kinks in the curve for A^(0, 7) vs. 7 no longer fall at the position of the nodes. 
This would be the case in a tetragonal system but is not so for the ellipsoidal Fermi surface 
of Fig. 0. We see from Eq. (^) that the kink occurs at a critical angle 7c given by (note 
that s = 0): 

1 -I- a 

tan7c = , (26) 

1 — a 

which is 7r/4 only for the case a = (no anisotropy in the electronic band structure). This 
is an important result because it is often believed that the sharp structure in A^(0, 7) as a 
function of 7 - the magnetic field orientation - gives the position of the nodes. This is not 
true when band anisotropy exists. 

The geometry of the situation is shown in Fig. ^ which depicts the laboratory frame while 
the result (p^ ) was obtained in the transformed, p-frame where the Fermi surface is a circle. 
In this figure k„ gives the position of a gap node on the ellipsoidal Fermi surface = 0, in 
Eq. (|l|). The Fermi velocity at that point, vpCkn), points in the direction tan~^ (^j^t&iKp^. 
If the magnetic field H is placed in that direction the circulating supercurrents will be in a 
plane perpendicular to the two dimensional Brillouin zone and oriented along the tangent 
to the Fermi surface at k = k„ as shown in the figure by the thick shaded line. Thus, there 
will be no Doppler shift. In the laboratory frame tanip = ^(1 + s)/(l — s) from Eq. (^ and 
we get immediately the result 

tan7,= - J- 27 

1 — a V 1 — s 

for a finite sub dominant s-wave part to the gap. This same result follows from Eq. (^4]) 
when the expression in the absolute value switches sign. For finite values of s, the angle 
7c coincides with the position of the nodes in the {d + s)-symmetric gap in the laboratory 
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frame only if a = 0. We can also see from (pTj) that a negative value of the subdominant 
s-wave part {i.e. a negative value of s) partially compensates for the effect of a but does 
not cancel it out completely. 

In Fig. ^ we show further results for the zero frequency QDOS A^(0, 7) as a function of 
magnetic field orientation 7 for increasingly anisotropic bands. The parameter a is increased 
from 0.2 in the sohd curve to 0.4 (dashed curve), 0.6 (dotted curve), and 0.8 (dashed- 
dotted curve). It is clear that the amplitude of the predicted oscillations in the QDOS 
increases considerably with increasing effective mass anisotropy (a) and that consequently 
this anisotropy will be easier to detect as the orthorhombicity of the band structure increases, 
i.e. as a increases. There is nearly a factor of 6 between the value at minimum and 
at maximum in the dashed-dotted curve. The pattern of oscillation is also changed as a 
increases. The intermediate secondary maximum, around 7 = 1.5 rad in the solid curve, is 
not present any longer in the dashed-dotted curve. Instead there is a minimum. 

A fit to data of the form (|2^) could in principle yield information on the two most 
important parameters a and s but the position of the sharp structure in the QDOS by itself 
does not give directly the position of the nodes or the ratio of the subdominant s-wave 
admixture to the dominant (i-wave component. It is the direction of H giving no Doppler 
shift which determines the position of the minimum in A^(0,7). Before leaving Eq. (^) we 
note that A^(0,7) is proportional to u = Eh/^o and this is the same dependence on the 
magnetic field as has been found in the tetragonal case. There is no explicit change in this 
quantity due to anisotropy so that the scaling variable for the QDOS and for the specific 
heat stays unchanged. This holds for the pure case but will be modified when disorder is 
considered. In that case scaling breaks downi and specific details play a role. 

The finite frequency QDOS in a pure orthorhombic superconductor with a (rf + s)- 
symmetric gap is also of interest. It follows from the formula 

^(.;)-vHp)q. 1 
v/(^(^)-v^(p)q,)2-A2^,(^,0)J 

which has been written in a form which remains valid even when impurity scattering is 

included. Here p = r/R. For the pure case uj{uj) = uj and K(Ij^s{uj , (p) = Ao/(0). Results of 
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the numerical evaluation of Eq. (^) (an approximate analytic formula similar to Eq. (p4D 
can also be derived but is not given here as it has the same form as already presented in 
the work of Vekhter et a/.lllll3) gives the results shown in Fig. ^ for four directions of the 
magnetic field. Here we present numerical results for N{uj,'-f) as a function of 00/ Aq for 
a = 0.2 and u = 0.3. For the magnetic field in the direction of 7 = 0° (antinodal direction 
in pure d-wave) the solid curve applies and an u'^ behavior is obtained at small u, just as 
in the tetragonal case. 00 But the linear dependence expected in the nodal direction in 
the tetragonal case now occurs in some other direction where the Doppler shift at the nodal 
position is zero, namely at 7c as defined by Eq. (^) which is where the superfiuid momentum 
qs is perpendicular to the Fermi velocity. For a = 0.2, 7c ^ 56° the dashed-dotted curve 
indeed displays linear dependence. 



III. IMPURITY EFFECTS 

We now include in the calculations elastic impurity scattering. In this case the full self 
consistent Green's function (H) enters the impurity term and we need to solve two coupled 
equations, one for the renormalized Matsubara frequency ijj{uj) and the other one for the 
gap Ad+si^, the real frequency axis they are: 

Lj{p, P,u) = uj + i7it^Q{p, P,u!), (29a) 
Arf+,(p, p, to, 0) = Ad+si(p) + int+D{p, p, to), (29b) 

with Ad+s{4') = Ap of Eq. (|^, the pure system gap in the transformed frame. Here p 
and (3 are associated with the vortex and is the polar angle defining the direction of 
p. The use of Eqs. (p9D confines our analysis to low temperatures where the temperature 
dependence of the gap function can be neglected. Moreover, these equations are only valid 
in Horn's approximation, or weak impurity scattering limit. A slightly more complicated set 
of equations needs to be solved in the unitary limit (strong impurity scattering), namely 

n^{p,(3,u;) + D^{p,p,u;) 
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Ad+s{p, P, uj, 0) = Ad+s{(p) + iirT' 



(30b) 



In the tetragonal case (a = 0), also considered here, the symmetry is described by the 
point group and the {d + is)-symmetry of the order parameter is the only possible mixed 
symmetry allowed under the symmetry operations of this group. In this case Eqs. (|29bD and 



( pOb|) are to be modified to give 



A,(p, P, Lj) = As + int+D{p, p, u) 



(31) 



in Born's limit and 



A,{p,p,Lj) = A, + inT^ 



D{p,P,u;) 



(32) 



n^ip,P,u;) + D^{p,P,u;y 

in the resonant scattering limit. is the s-wave contribution to the order parameter in the 
clean limit and is equal to zero in pure li-wave symmetry. The renormalized gap function is 
then written as: 



Ad+isip, P, UJ, 0) = Ad(0) + iA,{p, p, uj). 
Finally, the functions D{p,P,u) and fl{p,P,uj) are defined as 

t^(p./3,tx;)-vf (p)qs \ orthorhombic, 

tetragonal, 
orthorhombic. 



(33) 



n{p,p,uj) 



D{p,P,uj) 



(ii(p,/3,i^)-VF(p)qa) -A2^^(p,/3,w,9i)/ 
iI'(p,/3,t^)-VF(p)qa 

^d+s{P,l^,i^,<t>) ^ 



(34a) 



(cD(p,/3,c^)-VF(p)q.)^-A^^^(p,/3,c^,, 



(34b) 



tetragonal. 



, V'(i,(p,/3,a;)-v_p{p)q.)^-Ai(p,/3,^)-A2{<;i) 

with (■ ■ ■)(^ the Fermi surface average in the p-frame. Finally, Eq. ( ^81) still applies for the 
QDOS and it is only at this later stage that an average over vortex variables is carried out. 
Eqs. ( [25| ) or (^) need to be solved self consistently for different values of p, P, and u. The 
two remaining parameters are t'^ or r+ which are related to the strength of impurity scat- 
tering in the normal state where nt^ = ttF^ = l/2rjmp and Ti^p is the impurity scattering 
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time in the normal state. It should also be noted at this point that the impurity scatter- 
ing potential in the laboratory frame is now anisotropic. The role of anisotropic impurity 
scattering in anisotropic superconductors has been studied by Harah and Nagi.0 Accord- 
ing to their results the influence of the impurity concentration on superconductivity varies 
with the anisotropy of the scattering potential, qualitative changes in, say the temperature 
dependence of the penetration depth or the energy dependence of the optical conductivity, 
cannot be expected. 

A first result with impurities is given in Fig. |^ for Born scattering with = Eh/^Q in 
a tetragonal system with a pure d-wave gap. Here, Aq = 24meV, u = 0.1, and what is 
plotted is the zero frequency value of the QDOS, A^(0,7), as a function of the orientation 
of the in-plane magnetic field H which makes an angle 7 with the x-axis. The solid line 
includes impurities and is compared with the dotted line for the pure case. We note that 
the minima are somewhat filled and rounded by the impurities as we might have expected 
and the maxima are slightly higher and broader. Impurities increase the QDOS in the low 
energy region. Close examination of the two curves also shows that they both have four-fold 
symmetry. This is to be expected, but it is pointed out that, even for a tetragonal system, 
the introduction of impurities introduces a finite s-wave component to the gap, which is now 
of {d + zs)-symmetry, in Eqs. (^) or ( |32D because it is the self consistent Green's function 
that enters the definition of (|3^ and this leads to a nonzero value of D{p,f3,uj). Our final 
result remains four-fold symmetric, i.e. it retains the full tetragonal symmetry. 

The effect of impurity scattering is much larger in the unitary limit. Results are shown 
in Fig. ^. Again the underlying system has tetragonal symmetry and a pure d-wave gap. 
The impurity scattering parameter is F"^ = Eh/^Q with a gap amplitude Aq = 24meV and 
u = Eh/ ^0 = 0.1 in Eqs. (|5(]|). The solid line represents our results for the value of the 
QDOS A^(0, 7) as a function of the magnetic field orientation 7 in the CUO2 plane. Notice 
the scale on the vertical axis and that the size of A^(0, 7) is roughly three times its value 
for the Born limit (Fig. |^. This changes quite a lot when the frequency is increased. The 
dashed curve applies to cj/Aq = 0.05, the dotted one to cu/Aq = 0.1, and the dashed-dotted 
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one to cu/Ao = 0.15. The percent anisotropy in the QDOS is reduced with increasing oj. 



Results for the orthorhombic case are shown in Figs. ^ and [Ty. In these figures v = 
Eh/^o = 0.1, the band anisotropy a = 0.41, and the s-wave admixture to the gap is 
characterized by s = —0.25 taken from a previous fit to experimental data in YBCO as 
discussed by Wu et a/.il and Schiirrer et a/.il Reasonable values for the parameters in 
our admittedly simplified band structure model, namely a = {itlx — ruy) / {m.j. + niy) and 
the gap anisotropy, namely s, can be set in comparison with penetration depth data on 
YBCO at optimum doping. The measured ratio of the zero temperature penetration depth 
depends on the anisotropic effective masses and sets the value of a = 0.41 which follows from 

Aa/Afe = 1600A/1030 A.il'ili§il Furthermore, in our model the low temperature slopes of 

the penetration depth are approximated by 

■(1 — a — s), 

■(1 + a + s). (35) 
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The experiments of Bonn et al£3 give s = —0.25 as an estimate of the orthorhombicity in 
the gap. It is these values that we have used to estimate the expected anisotropy in the 
QDOS for H in the Cu02 plane. As we have already stressed in the introduction, the band 
structure in YBCO is much more complex than the model treated here but the simplicity of 
our model allows us to get insight into the role of orthorhombicity and develop a qualitative 
picture of its effect. There are no adjustable parameters left. 

For the orthorhombic case the pattern in Fig. ^ for A^(0, 7) vs. 7 is two-fold symmetric 
even in the pure case. As before. Born scattering (solid curve with t+ = Eh/1Q) simply 
smooths out slightly the minima seen in the pure case curve (dotted line). On the other 
hand, unitary scattering, solid gray curve for r+ = Eh/ 16, pushes the curve up by a factor 
of about two over the pure limit case and the pattern is also significantly modified and 
smoother. The small maxima around 7 = 1.5 rad and symmetry related points in the pure 
curve (dotted) are completely washed out by unitary impurity scattering. The amplitude of 
the variations in value of A^(0, 7) is also considerably reduced. 

The curves described so far apply to untwinned single crystals. Three more curves are 
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shown in Fig. |^. The dashed-dotted curve apphes to the case of twinned samples and was 
obtained from the clean limit (dotted curve) by averaging it with a similar curve displaced 
by 90°. This corresponds to a simple average of equal numbers of twins with x- and y- 
axis interchanged. It is clear that this average greatly reduces the predicted anisotropy and 
experiments aimed at discovering this anisotropy are best made on untwinned samples. This 
is confirmed in preliminary data from Junod's groupil who indeed find no significant in-plane 
anisotropy in the specific heat on twinned single crystals. At first sight the noise in their 
experiment is below, but of the order of the anisotropy predicted in our calculations. There 
are many reasons why the actual anisotropy might in fact be even less than the amount 
shown in our pure twinned curve (dashed-dotted). To mention only one: YBCO is more 
three dimensional than many of the high oxides and this implies one more integration in 
the definition of the QDOS (|2^) and this can be expected to reduce the anisotropy further. 
The two other curves in Fig. ^ include impurities and show even less anisotropy than does 
the pure case. The dashed curve applies to Born scattering with t+ = Eh/^Q and the dashed 
gray curve to resonant scattering with F"^ = Eh/^Q- 

In Fig. IIU] we present finite frequency results for N{uj,'-)) vs. 7 (the magnetic field 
orientation in the CUO2 plane) for four finite frequency values, namely cu/Aq = 0.05 (dashed- 
dotted curve), tu/Ao = 0.1 (dashed curve), tu/Ao = 0.15 (dash-double dotted curve), and 
cij/Ao = 0.2 (dotted curve). The results are to be compared with the solid curve which 
applies to the case = and has already been shown in Fig. |^. It is clear that the two 
pronounced minima in the solid curve become shallower with increasing uj and are already 
quite small for the dotted curve with cj/Aq = 0.2. This complicated pattern should be 
reflected in the accompanying pattern for the temperature dependence of the specific heat. 

More insight into the frequency dependence of the anisotropy can be obtained from the 
frequency dependence N{uj, 7) which is shown for several fixed directions 7 in Fig. |11], frame 
(a) and (b). These apply respectively to the tetragonal and orthorhombic case. In frame (a) 
7 = 0° (antinodal direction) and 7 = 45° (nodal direction) are considered. For the pure case 
the solid curve which applies for 7 = 0°, shows an u;^ dependence at small uj while the dotted 
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curve for 7 = 45° is linear as expectedtJo and falls below the value along the antinodal 
direction. The nodal direction is indeed the direction of the minimum in A^(0,7) of Fig. |^. 
As the frequency uj is increased the solid and dotted curves rapidly come together and even 
cross so that the nodal direction has a larger value of A^(u;,7) than the antinodal direction 
for uj larger than the crossover frequency. The anisotropy is now very small and will not be 
observable in specific heat experiments. Dependencies of A^(co',7) vs. uj are changed only 
slightly with the introduction of Born scattering impurities; the dotted curve is for 7 = 0° 
and the dashed-dotted one for 7 = 45°, and t"*" = Eh/^Q- On the other hand, the gray curves 
for resonant scattering with F"*" = Eh/^Q are very different and are approximately constant 
at low UJ with little difference between the two directions. It is obvious that little anisotropy 
remains in the QDOS. Similar results are given in frame (b) for the orthorhombic case with 
a = 0.41, s = —0.25, and u = 0.1, i.e. the same value of Eh to gap amplitude is used as 
in frame (a). The most important difference to be noted is that the critical angle, at which 
an u dependence is found in the clean limit, is at the angle 7 = 61.3° and not at 7 = 45°. 
Also, the anisotropy is not washed out quite as much in the resonant scattering case when 
compared to the tetragonal case. Further in this case at sufficiently high frequencies the 
two curves for 7 = 0° and 7 = 61.3° cross and the curve for 7 = 0° falls below the one 
for 7 = 61.3° thus reversing the pattern of maxima and minima shown in Fig. ^ for lower 
values of u. We also note that as u increases out of zero the QDOS N{uj, 7) decreases for 
both 7 values given in frame (a) (tetragonal case) while for the orthorhombic case in frame 
(b) the reverse is true. This effect is also seen in Figs. ^ and In the tetragonal case 
the solid line of Fig. ^ for which = is above the others for u ^ while in Fig. 0, the 
orthorhombic case, it is below. 



IV. CONCLUSIONS 

We have calculated the effect of the circulating vortex supercurrents on the QDOS of a 
two dimensional d-wave superconductor. For the magnetic field H placed in the CUO2 plane 
the QDOS varies as a function of the angle 7 between H and the x-axis. We have included in 
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our calculations an a-b effective mass anisotropy of about 2 with a view of modeling the band 
structure of YBCO which is orthorhombic because of the existence of CuO chains, and for 
which the charge carrier oscillator strength on the chains is approximately of the same order 
as in the Cu02 planes. In addition, the gap on the Fermi surface was assumed to posses a 
subdominant s-wave as well as a dominant d-wave character. Consideration of penetration 
depth data fixes the value of this admixture so that we have no adjustable parameters. 
Impurities are also included in our calculations in both, Born's and unitary limits. It is 
found that the four-fold pattern predicted for the dependence of the zero frequency value 
of the QDOS in a pure tetragonal system is now changed to a two-fold variation due to the 
orthorhombicity in the gap or in the band structure. The introduction of impurities leaves 
the symmetry of the pattern unchanged but can reduce its amplitude. It is also found that 
the pattern of variation of N{u!, 7) with 7 is altered as the frequency is increased, although 
it remains four-fold for tetragonal systems and two-fold for the orthorhombic case. 

An important result is that the angle at which minima occur in iV(0, 7) vs. 7 does 
not correspond to the angle at which the zeros occur in the gap as would be the case in a 
tetragonal system. Instead, the minima occur when the normal state quasiparticle Fermi 
velocity in the nodal direction is perpendicular to the direction of the vortex supercurrents, 
i.e. parallel to the direction of the external magnetic field. 

Specific heat measurements in clean samples are most favorable for observing the pre- 
dicted two-fold pattern of anisotropy. To observe it in YBCO detwinned samples are best 
because an average over equal numbers of twins greatly reduces the predicted anisotropy for 
detwinned samples. Consideration of more realistic tight binding band structure models are 
in progress. Our very simple effective mass model, however, has allowed us to understand 
the important changes orthorhombicity brings to the magnetic field orientation dependence 
of the QDOS for a two dimensional CUO2 plane with H in the plane. It has helped us to 
determine the experimental conditions most favorable to the observation of the predicted 
anisotropics. 
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FIGURES 

FIG. 1. The ellipsoidal Fermi surface in k-space given by ek = in Eq. (|^) (solid line). Here we 
have used a = 0.4. Also shown is the vector k with the direction ip and the corresponding Fermi 
velocity at k, Vi;'(k). Note that it is not parallel to k. The transformed Fermi surface to p-space 
is a circle and is shown as a solid, gray line, pp = y/2rh£F, with fh = (nix + my)/2 and the angle 
(p is identified. 

FIG. 2. Order parameters on the transformed Fermi surface (p-frame) as given by Eq. (^). (a) 
a = s = (pure d-wave gap); (b) a = 0.4, s = 0; (c) a = s = 0.4; (d) a = 0.4, s = —0.4. The 
straight lines denote the nodal angles. 

FIG. 3. The zero frequency value of QDOS, A^(0, 7), as a function of magnetic field orientation 
in the Cu02 plane with 7 the angle of H with respect to the kx-axis. The solid line is based on the 
approximate analytic formula (|2^ while the dashed line is obtained from numerical evaluation of 
the exact one dimensional integral ([l8|). Only band anisotropy (m^. it niy) is included with a = 0.2 
(ellipsoidal Fermi surface) and the gap is pure d-wave with nodes in the laboratory, (k)-frame at 
ib7r/4 and symmetry related points. In the figure u = Eh/^o = 0.3. 

FIG. 4. The Fermi surface (ellipse) in the laboratory frame. Also shown is the direction of the 
node k„ of the gap {(fn)- The direction of the magnetic field H with the angle 7 is taken to be 
parallel to the Fermi velocity at the node Vir(k„) so that the supercurrents defining the vortex 
flow in a plane (indicated by the thick, solid, gray line) are perpendicular to the kx-ky plane and 
parallel to the tangent on the Fermi surface at Vi?(k„). 

FIG. 5. The QDOS A^(0, 7) at zero frequency as a function of magnetic field orientation in 
the Cu02 plane (7) for various values of the band structure anisotropy (a). As a increases the 
amplitude of the oscillations increases: the solid line is for a = 0.2, the dashed one for a = 0.4, 
the dotted one for a = 0.6, and, finally, the dashed-dotted one for a = 0.8. 
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FIG. 6. The frequency dependence (a;) in units of Aq (w/Aq) of the QDOS, N{uj,^) for the 
pure case at four different magnetic field orientations, namely 7 = 0° (solid line), 45° (dotted line), 

51° (dashed line), and 56° (dashed-dotted). The 7 = 0° curve shows behavior at low uj while 
the 7 = 56° is linear in uj as would be the case in zero field for a d-wave superconductor. 

FIG. 7. The zero frequency value of the QDOS A'^(0,7) as a function of the magnetic field 
orientation 7 for a pure d-wave tetragonal {a = s = 0) superconductor. The solid line includes 
some impurities characterized by = Eh/^G and the dashed curve is the clean limit result. 

FIG. 8. The dependence of the QDOS, N{aj,^) on the orientation of the magnetic field in the 
Cu02 plane (7) for four different frequencies u) in the unitary limit with = Eh/^Q- The solid 
curve is for w = 0, the dashed one for w/Aq = 0.05, dotted for w/Aq = 0.1 and dashed-dotted for 
a;/Ao = 0.15. u = Eh/^o = 0.1 in a tetragonal system with a pure d-wave gap. 

FIG. 9. The dependence of the zero frequency QDOS, N{0, 7), on the magnetic field orientation 
7 in the CUO2 plane for an orthorhombic system with a = 0.4, s = —0.25, and u = Eh/^o = 0.1. 
The dotted curve is for the pure system, the solid one for Born scattering with i+ = Eh/IQ, and 
the solid gray curve is for unitary scattering with r+ = Eh/IG- The dashed-dotted curve was 
obtained from the dotted one by taking the average of it, as is, and its displacement by 90°; it 
applies to experiments that average over equal numbers of twins. The dashed curve includes Born 
scattering and the dashed gray curve unitary scattering. Both are for twinned crystals. 

FIG. 10. The dependence of the QDOS, N{uj,'~f), on the orientation of the magnetic field 
7 in the CUO2 plane for different frequencies, namely uj/Deltao = (solid line), ui/Aq = 0.05 
(dashed-dotted hne), lo/Aq = 0.1 (dashed line), iv/Aq = 0.15 (dashed-double dotted line), and 
lo/Aq = 0.2 (dotted line). Orthorhombic d-wave with a = 0.4, s = —0.25, v = 0.1, and 

r+ = Eh/iq. 



23 



FIG. 11. The frequency dependence {u) in units of the gap amphtude Aq of the QDOS, N{u, 7), 
for various values of the magnetic field orientation 7 and impurity content. The solid and dashed 
lines arc for the pure case. The dashed and dashcd-dotted ones are for = E^/^Q and the solid 
gray and dashed gray ones for r+ = Eh/1Q- In each pair of curves the first is for 7 = 0° while the 
second is in the direction where two node regions have no Doppler shift due to the magnetic field. 
In (a) which applies to the tetragonal case this is for 7 = 45° which is the nodal direction but in 
(b) which applies to the orthorhombic case the critical direction is for 7 = 61.3°. This is different 
from the nodal direction in the laboratory frame. 
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